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ABSTRACT 

Using high resolution cosmological simulations, we study hydrogen and helium gravitational cooling 
radiation. We focus on the He II cooling lines, which arise from gas with a different temperature 
history (T max ~ 10 5 K) than H I line emitting gas. We examine whether three major atomic cooling 
lines, H I A1216, He II A1640 and He II A304, are observable, finding that Lya and He II A1640 cooling 
emission at z — 2 — 3 are potentially detectable with deep narrow band (R > 100) imaging and/or 
spectroscopy from the ground. While the expected strength of H I A1216 cooling emission depends 
strongly on the treatment of the self-shielded phase of the IGM in the simulations, our predictions 
for the He II A1640 line are more robust because the He II emissivity is negligible below T ~ 10 4 5 K 
and less sensitive to the UV background. Although He II A1640 cooling emission is fainter than Lya 
by at least a factor of 10 and, unlike Lya, might not be resolved spatially with current observational 
facilities, it is more suitable to study gas accretion in the galaxy formation process because it is 
optically thin and less contaminated by the recombination lines from star-forming galaxies. The He II 
A1640 line can be used to distinguish among mechanisms for powering the so-called "Lya blobs" - 
including gravitational cooling radiation, photoionization by stellar populations, and starburst-driven 
superwinds — because (1) He II A1640 emission is limited to very low metallicity (\og(Z/Z Q ) < —5.3) 
and Population III stars, and (2) the blob's kinematics are probed unambiguously through the He II 
line width, which, for cooling radiation, is narrower (a < 400 km s _1 ) than typical wind speeds. 
Subject headings: cosmology: theory — galaxies: formation — intergalactic medium 
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1. INTRODUCTION 

Galaxies grow partly by accretion of gas from the sur- 
rounding intergalactic medium and partly by mergers 
with other galaxies. Observational studies of galaxy as- 
sembly have focused primarily on merger rates, which 
can be measured indirectly by counting close pairs and 
merger remnants. However, all the mass that enters 
the galaxy population ultimately does so by accretion 
— mergers can only redistribute this mass from smaller 
systems to larger systems. Furthermore, numerical sim- 
ulations predict that even large galaxies grow primar- 
ily by smooth gas accr etion rather than by cannibal- 
ism o f smaller objects l|Murali et al.l 120021 iKeres et al.l 
12005^1 . Gas shock- heated to the virial temperature of a 
typical dark matter halo would radiate most of its ac- 
quired gravitational energy in the soft X-ray continuum, 
making individual sources very difficult to detect, espe- 
cially at high redshift. However, Fardal et al. (2001, 
hereafter F01) show that much of the gas that enters 
galaxies in hydrodynamic cosmological simulations never 
heats to high temperatures at all, and that it therefore 
channels a substantial fraction of its cooling radiation 
into atomic emission lines, especially H I Lya. F01 and 



lHaiman. Soaans. fc Quataertl (120001) sugge sted that ex- 
tended "Lya blobs" (e.g.jKeel et al]|1999t ISteidel et al 



1 20001 iFrancis et all200^lMatsuda et al.l2004UDev et al 

2005), with typical sizes of 10 — 20" and line luminosi- 
ties LLya ~ 10 44 ergs s _1 , might be signatures of cooling 
radiation from forming galaxies. iFurlanetto et al.l 1)20051) 
have also investigated predictions for Lya cooling radia- 
tion from forming galaxies in hydrodynamic simulations. 

In this paper, we investigate other aspects of cooling 
radiation from forming galaxies, in particular the po- 
tentially detectable radiation in the He II A304 ( "Lya" ) 
and He II A1640 ("Ha") lines of singly ionized helium. 
While challenging, the successful detection He II line 
emission would complement H I Lya measurements in 
at least three ways. First, because H I and He II 
line cooling rates peak at different temperatures (T ~ 
10 4 3 K vs. T ~ 10 5 K), measurements of both lines 
could constrain the physical conditions of the emitting 
gas. Recent theoretical studies imply that "cold mode" 
accretion, in which the maximum gas temperature is 
well below the halo virial temperature, is a ubiqui- 
tous and fundame ntal feature of galaxy formation (F01; 
iKatz et all 12003 iBirnboim fc Dekell l200l IKeres et alJ 
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20051 iDekel fc Bimboiml200l. a v iew an ticipated by the 
early analytical work of lBinnevI l|1977fl . Simultaneous 
measurements of He II and H I emission from different 
types of galaxies could eventually test detailed predic- 
tions for the temperature distribution of cooling gas. 

The second advantage of the He II A 1640 line is that 
it should be optically thin, allowing a straightforward 
interpretation of the observed emission in terms of the 
spatial distribution and kinematics of the cooling gas. In 
contrast, the radiative transfer effects on the H I Lya 
emission from accreting gas are more complicated, an 
issue that we will investigate i n future work (J. Kollm eier 
et al., in preparation; see also lCantalupo et al.ll2005lL 

Third, the different temperature dependence and ab- 
sence of radiative transfer effects in He II lines could help 
distinguish cooling radiation from alternative explana- 
tions of "Lya blobs," such as emissi on from collisionally 
ionize d gas in galactic superwinds (jTaniguchi & Shiova 
2000) or from gas photoionized by young stellar popu- 
lations. For example, only the lowest metallicity stars 
(logZ/Z© < —5.3) have hard enough spectra to ionize 
He II to He III , so stellar photoionization will generally 
not produce He II line emission. 

It is possible that He II cooling emission from forming 
galaxies at high redshift has already been detected. One 
such example is the broad He II A1640 line in the com- 
posite spectrum of L yman break galaxies (LBGs) from 
iShaplev et al.l l)2003|) . The composite He II line shows a 
rather broad line width (FWHM~ 1500 km s -1 ), which 
is a possible signature of Wolf-Rayet stars. However, it 
is difficult to reproduce the strength of the He II line 
via s tellar population m odels with reasonable parame- 
ters l|ShaDlev et al]l2003lL We show in this paper that 
the He II cooling emission around individual galaxies is 
detectable, which suggests it might be fruitful to search 
for He II A1640 cooling emission in the outer parts of 
individual LBG's. 

The next section describes our simulations and the ra- 
diative cooling mechanisms that lead to H I A1216, He II 
A1640, and He II A304 line emission. We present H I 
A1216 and He II A1640 cooling maps in t !3.1l and dis cuss 
the properties of the cooling radiation sources in H3.2I 
In H3.3I we examine the detectability of the three major 
cooling lines in the far ultraviolet and optical and discuss 
the best observational strategies. In 21 we discuss mech- 
anisms other than gravitational cooling radiation that 
can generate extended Lya emission and describe how 
the He II cooling line might help us discriminate among 
those mechanisms. We summarize our conclusions in SjSJ 

2. SIMULATIONS AND COOLING RADIATION 

We use Parallel TreeSP H simulations 
ijDave. Dubinski. fc Hernauistl 11997(1 including the 
effects of radiative cooling, star formation, thermal 
feedback, and a spatially uniform metagalactic pho- 
toionizing background. We analyze two simulations: 
one with a cubic volume of 11.111 h^ 1 Mpc (comoving) 
on a side and a spatial resolution of 1.75 hr 1 kpc 
(comoving; equivalent Plummer softening), the other 
with a cubic volume of 22.222 h" 1 Mpc on a side 
and 3.5 h^ 1 kpc resolution. Hereafter, we refer to 
these two simulations as the 11 Mpc and 22 Mpc 
simulations, respectively. The simulations consist of 
128 3 dark matter particles and 128 3 gas particles, 



giving a mass resolution of jtisph = 1-3 x 10 7 M Q 
and mdark = 10 8 Af© for the 11 Mpc simulation, and 
tosph = 1-1 x 10 8 M Q and m dark = 7.9 x 10 8 Af Q for the 
22 Mpc simulation. We adopt a ACDM cosmology with 
the parameters J7m — 0.4, J7a = 0.6, fif,/i 2 = 0.02, and 
Hq = 65 km s _1 Mpc -1 . In our calculation of cooling 
emission lines, we basicall y follow the radiative cooling, 
proces ses described in iKatz. Weinberg, fc Hernauistl 
(1996, hereafter KWH). Here, we briefly summarize the 
cooling processes that can contribute to line emission. 

The four underlying assumptions of radiative cooling 
are: primordial composition, ionization equilibrium, an 
optically thin gas, and a spatially uniform radiation field. 
In the simulations, we adopt X — 0.76 and Y = 0.24, 
where X and Y are the hydrogen and helium abundance 
by mass. The abundances of each ionic species (H°, H + , 
He , He + , He ++ ) are solely determined by assuming that 
the primordial plasma is optically thin and in ionization 
equilibrium (but not in thermal equilibrium). The func- 
tional forms of the temperature-dependent recombina- 
tion rates, collisional ionization rates, collisional excita- 
tion rates, and the rate equations are given in §3 of KWH 
and the tables therein. The uniform photq i onizin g UV 
background is taken from lHaardt fc Madaul l|1996[) . 

In the next section, we describe the Lya and He II 
cooling curves under the optically thin gas assumption, 
the assumption used in our simulations. However, in the 
high density regime, a gas cloud becomes dense enough 
to shield the central part of itself from the UV back- 
ground, i.e., becomes self-shielded. Therefore, the actual 
emissivity of this self-shielded (or condensed phase) gas 
is highly uncertain. We discuss below how to treat this 
phase of the IGM to derive better estimates of the cooling 
radiation. Once we generate the Lya and He II cooling 
emissivities, the cooling radiation is determined by how 
many gas particles populate a certain range of tempera- 
ture and density, where we appeal to the high-resolution 
cosmological simulations mentioned above. 

2.1. Cooling Curves 

Among the various radiative cooling processes, only 
two can produce H I Lya (1216A), He II Lya (340A), 
and He II Balmer a (1640A) photons: the recombination 
cascades of a free electron and the collisional excitation of 
a bound electron to an excited state followed by radiative 
decay. The dominant cooling mechanism is the collisional 
excitation of neutral hydrogen and singly ionized helium, 
which have their peaks at temperatures of T ~ 10 4 ' 3 K 
and ~ 10 5 K, respectively. Figure^shows the H I A1216 
and He II A304 emissivities for gas of different densities in 
ionization equilibrium in the presence of a metagalactic 
photoionizing background. The dashed lines represent 
the collisional excitation cooling lines of neutral hydro- 
gen and singly ionized helium. The dot-dashed lines and 
dotted lines denote the recombination lines for these two 
species due to photoionization and collisional ionization, 
respectively. The solid lines represent the total Lya cool- 
ing rates of hydrogen and helium. Below T ~ 10 4 K, 
collisions with free electrons are not energetic enough to 
raise bound electrons to upper levels or to ionize the neu- 
tral hydrogen, so the collisional cooling rate of hydrogen 
drops quickly below this temperature. For singly ionized 
helium, the collisional cooling rate drops to virtually zero 
below T ~ 10 4 - 6 K. Therefore, below T ~ 10 4 K, pho- 
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Fig. 1. — Normalized line emissivity logA/n^ as a function of temperature for a primordial plasma at densities p/pb = 10 2 , 10 3 , and 
10 4 (from left to right) in the presence of a UV ionizing background at z = 3. In each panel, the dashed lines represent the collisional 
excitation cooling rates. The dot-dashed lines and dotted lines denote the recombination rates owing to photoionization and the collisional 
ionization, respectively. The solid lines represent the total line cooling rates of hydrogen and helium. The bold solid lines below the He II 
cooling curves represent the He II A1640 line emissivity. Compared with H I, the cooling rates of He II owing to the UV ionizing background 
become significantly weaker as the gas density increases. 



toionization by the background UV spectrum (heating) 
and the following recombination (cooling) is the domi- 
nant source for cooling line emission of the primordial 
plasma. 

In the presence of photoionization, the cooling curve 
depends on the density as shown in Figure ^ because 
the equilibrium abundances of each species depend on 
the density. While the collisional ionization rate per 
unit volume scales as n^, the photoionization rate per 
volume is proportional to only nn- In low density gas 
(p/pb IS 10 2 ), a significant fraction of H I and He II 
is photoionized and their abundances are mainly deter- 
mined by the photoionization process, so the collisional 
excitation feature is relatively weak. As the gas density 
increases (p/p~b ~ 10 3 — 10 4 ), however, collisional excita- 
tion becomes more important and the cooling curves ap- 
proach pure collisional equilibrium, the so-called coronal 
equilibrium. The notable difference between the cooling 
curves of hydrogen and singly ionized helium is that the 
cooling rates of He II owing to photoionization become 
significantly weaker as the gas density increases. If we 
assume that the medium is optically thin to the ioniz- 
ing background — i.e., not self-shielded — hydrogen is 
almost fully ionized over the entire temperature range 
even at high densities, but ~90% of the helium is in a 
singly ionized state below T ~ 10 4 8 K at high densities. 
Thus the assumption that the gas is optically thin ev- 
erywhere is roughly valid for the calculation of the He II 
line fluxes, while it is poor for hydrogen. The correction 
for the self-shielding of hydrogen will be discussed in the 
next section. 

We estimate the He II A1640 flux from the He II A304 
flux by considering the ratio of He II A1640 to He II A304 
in the recombination cascades and in the collisional ex- 
citations, respectively. The thick solid lines below the 
He II cooling curves in Figure ^ represent our estimate 
of the He II A1640 line emissivity. Below T ~ 10 5 K, the 
optical depth of He II A304 is so large that case-B recom- 
bination is a good approximation. Even though the Lya 



optical depth is extremely large, the population of the 2p 
and 2s states is always much smaller than that of the Is 
state, because the de-excitation time for level transitions 
is very short (A2 P i s ~ 10 10 s _1 ). One might be concerned 
whether the population of the 2s state is large because of 
the forbidden transition (2s — > Is), but the two photon 
decay process is fast enough to de-populate 2s electrons 
(^2sis — 8.22Z 6 s _1 ). Therefore, Balmer lines are always 
optically thin. We adopt FJ^/F^a ~ 10% by extra po- 
lating the case-B values of IStorev &: Hummer! fl99 5) to 
the low density limit. For collisional excitation, we esti- 
mate the He II A 1640 flux using 

-^1640 ^ C±s3s + C\sZp + C\sZd hv\§AQ .j. 
•^304 C\ s 2p + Cis3 s + C\s3d ^304 

where Cy is the collisional excitation ra te from the i to 
the j state. We adopt the dj's from Ag garwal et alJ 
GUIEI- F i64p/ F 304 is roughly 2 - 4% in the temperature 
range of 10 5 < T < 10 5 7 where He II A304 collisional 
excitation cooling is dominant. In summary, the He II 
A 1640 flux is calculated by 

^Hcll A1640 = 

/rccn e nHc++aHc+^Lya + /coll"-e«Hc+ Cl2 hlS Lya , (2) 

where aS c + is the recombination rate into n ^ 2 states 
of He II, and by assuming / rcc w 10% and / co n - 4%. 

2.2. Self- shielding Correction 

A major di fferen ce between our work and that of 
iFardal et all (|2001^l is that our simulations include 
a uniform UV backg round radiation field (see also 
iFurlanetto et al.l[2005jT . However, because even state-of- 
art cosmological simulations like ours do not include ra- 
diative transfer, the self-shielded phase of the gas at high 
column densities is not treated properly. When the gas is 
heated to high temperature (T ~ 10 5 — 10 6 K) by falling 
into the forming galaxy's halo, the gas is mostly ionized, 
so it is reasonable to assume that the gas is optically thin 
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to the uniform UV background. Subsequently, when the 
gas cloud starts losing its thermal energy via cooling ra- 
diation, its neutral column density becomes sufficiently 
high that the metagalactic UV radiation cannot pene- 
trate the surrounding gas, and the cloud becomes self- 
shielded. 

Once the supply of ionizing photons is shut off, what 
happens to the self-shielded high column density clouds? 
First, the ionization states will achieve collisional ion- 
ization equilibrium, where the emissivity is determined 
solely by collisional ionization and collisional excitation. 
Second, because the cooling emissivity is boosted by 
these processes, the self-shielded cloud will cool more 
rapidly to T ~ 10 4 K than in the presence of heating 
by ionizing photons. Below this temperature, the subse- 
quent cooling is dominated by metal lines (if there are 
metals). Stars ultimately form from this cold gas. 

Because our simulations do not include the time evo- 
lution of the self-shielded gas or metal-line cooling, it 
is not clear how long the gas particles stay in the self- 
shielded phase and emit in collisional ionization equilib- 
rium. The Lya emissivities shown in Figure ^ become 
unreliable in this self-shielded regime. Thus we apply a 
pseudo self-shielding correction to the high density gas 
particles to correct their emissivities. This correction 
is not rigorous; to properly calculate the emissivity of 
the self-shielded phase of the IGM, one should incorpo- 
rate a radiative transfer calculation that includes non- 
uniform and anisotropic UV radiation fields. A different 
prescription for th e self-shielded phase is d efinitely pos- 
sible. For example, Furlanct to et alJ l)2005f) consider two 
extreme cases: 1) adopting zero emissivity and 2) using 
the collisional ionization equilibrium emissivity for the 
self-shielded phase. Our self-shielding correction scenario 
described below lies between these two extremes. 

To apply the self-shielding correction, we first define 
the "local" optical depth for each gas particle, T\ oca \(v) — 
Y^i n i (T i{ l/ ) (al), where n, and Oi are the number densi- 
ties and the photoionization cross sections of each species 
(H°, He , He + ), respectively. The "local" size of the gas 
cloud I — the length that corresponds to the volume 
that the gas particle would occupy in space — is de- 
fined as (mgas/p) 1 / 3 . 1 For each gas particle, the UV 
background spectrum J(v) is attenuated using this lo- 
cal optical depth, i.e. J(y)e~ Ttyl, \ and new photoioniza- 
tion/heating rates and equilibrium number densities are 
calculated. We then determine a new Ti oca i(^) from these 
values and iterate this procedure until the photoioniza- 
tion rates and the optical depths converge. We use these 
final ionization/heating rates to calculate the Lya and 
He II emissivity of each gas particle. This modified emis- 
sivity for each gas particle is what we will refer to as the 
self- shielding correction case. 

Clumping inside a gas particle and/or among gas particles 
could be approximated using a free parameter a such that al rep- 
resents the effective geometrical edge-to-center distance of the gas 
cloud. For example, a sp here = (3/4-7T) 1 / 3 is given for a single 
spherical gas cloud, whereas a = 2 corresponds to the clumping 
of (a/a sp hcrc) 3 — 34 gas particles. The value of a should vary 
from one particle to another, but we adopt a = 1 throughout the 
paper as a fiducial value. The over-density where the self-shielding 
occurs also depends on the choice of the edge-to-center distance al 
(and on the redshift). However, we find that the effects of adopting 
a = 0.5 — 2.0 is insignificant. 




log (p/pj 

Fig. 2. — The local optical depth of each gas particle as a 
function of over-density at z = 2. Solid and open circles indicate 
the optical depth at the H I (13.6 eV) and He II (54.4 eV) ioniza- 
tion edges, respectively. The horizontal and vertical lines represent 
t — 1 and p/pb — 10 , respectively. Note that H I optical depth 
increases abruptly at p/pb — 10 3 and the gas becomes self-shielded. 



In Figure |3 we show the local optical depths of each 
gas particle in the final equilibrium state as a function 
of over-density. For the 22 Mpc simulation at z = 2, we 
show t(H I) and r(He II), the optical depth at the H I 
(13.6 eV) and He II (54.4 eV) ionization edges, respec- 
tively. As indicated by the dotted lines, the H I optical 
depth increases abruptly from r ~ 1 to r = 10—100 
at p/pb — 10 3 . Therefore, the optically thin UV back- 
ground assumption is valid below p/pb = 10 3 , but the 
gas becomes self-shielded quickly above this over-density. 
Because the transition from the optically thin case to 
the self-shielded phase occurs abruptly, we also consider 
the emissivity of a condensed phase cut case as the most 
conservative for the cooling radiation. There we set the 
emissivity of the self-shielded gas particles to zero. 

Hence, in the following analyses, we consider three 
possibilities. First is the optically thin case that as- 
sumes a spatially uniform UV background for every gas 
particle. Second is the self- shielding corrected case de- 
scribed above that uses an attenuated UV background for 
each gas particle appropriate for the local optical depth. 
Third is the condensed phase cut case where we set the 
emissivity of gas with log T < 4.5 and p/pb > 10 3 to 
zero. We emphasize again that while none of these possi- 
bilities are rigorously correct, they range from the most 
optimistic (1) to the most conservative case (3). Note 
again that case (1) is appropriate for He II, but the full 
range of cases should be considered for H I. 

3. RESULTS 

3.1. Cooling Maps 

We generate HI A1216 and He II A 1640 cooling maps 
at z = 2 and 3 by applying our line emissivities to each 
pixel element and integrating them along the line of sight. 
The temperature and density of each volume element is 
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Fig. 3. — H I A1216 (top) and He II A1640 (bottom) cooling maps for the 11 Mpc simulation at Z = 3. The line of sight depth is 
Az ~ 0.019. The left panels show a part (i X 1) of our simulation, the middle panels show the brightest region at a finer pixel scale 
(~ 1.5 h _1 kpc per pixel), and the right panels show the cooling maps convolved with a 0"5 FWHM Gaussian filter to mimic a typical 
ground-based observation (re-binned to 0'.'2 per pixel). Note that we include Lya and He II emission from the IGM assuming that the 
emissivity of the condensed phase is zero, the most conservative case. The Lya cooling radiation from the gas around the forming galaxies 
will be observed as a diffuse and extended blob above ~ 10 -18 ergs s _1 cm -2 arcsec -2 , the current flux limit of ground-based detections 
(R = 100), whereas He II will be almost point source-like at current detection limits. 



computed using the usual SPH smoothing kernels, and 
the abundances of ionic species are calculated from these 
smoothed quantities. The thickness of the 11 Mpc sim- 
ulation along the line of sight is Az ~ 0.013 and 0.019 
for z = 2 and 3, respectively. We convert these cooling 
maps into surface brightness maps using our adopted cos- 
mology. Figure shows the H I A1216 and He II A 1640 
cooling maps for the 11 Mpc simulation at z = 3. We 
show the cooling maps for the condensed phase cut case 
to represent the most conservative prediction. The left 
panels show a part Q x j) of our simulation where the 
filamentary structure of the IGM — the so-called "cos- 
mic web" — is evident. In the middle panels, we show 
the brightest region (also the most over-dense region for 
z = 3) at a finer pixel scale (~ 1.5 hr 1 kpc per pixel; 
half of the spatial resolution of our simulation). The 
right panels show the cooling maps convolved with a 0'/5 
FWHM Gaussian filter to mimic a typical ground-based 
observation. 

As shown in Figure the Lya cooling emission from 
the IGM is somewhat extended above a surface bright- 
ness threshold of 10~ 18 ergs s _1 cm -2 arcsec -2 (the cur- 



rent limit of ground-based detections), whereas the He II 
emission will be seen almost as a point source. We will 
refer to these extended Lya cooling sources in our sim- 
ulations as Lya blobs hereafter. In the cooling maps, 
we include only Lya emission from the IGM, not from 
star formation (i.e. from photoionization caused by mas- 
sive stars followed by recombination). However, we find 
a compact group of stars and/or star-forming particles, 
i.e., galaxies, at the center of each Lya blob. There- 
fore, what we would actually observe are galaxies (or 
Lya emitters if dust absorption is negligible) embedded 
within the Lya blobs. 

The most uncertain factor in generating the cooling 
maps is how much the self-shielded gas contributes to 
the emission. To quantify this factor globally in the sim- 
ulations, we consider the H I A1216 and He II A 1640 
luminosity-weighted temperature and density plots for 
the optically thin (optimistic) case in Figure 0] To make 
these phase diagrams, we extract temperature and den- 
sity profiles for 100x100 evenly-spaced lines of sight, 
apply our emissivities to each radial bin, and integrate 
the temperatures and densities with the H I A1216 and 



6 



Yang et al. 



He II A1640 



6 - 



o 




H I A1216 




log (p/Pb) 



log (p/Pb) 



Fig. 4. — "Luminosity-weighted" phase diagram for 100 X 100 lines of sight through the simulation at z = 3 for He II A1640 (left) 
and H I A1216 (right). The diamonds indicate the lines of sight that have He II A1640 or H I A1216 surface brightnesses larger than 
10 — 19 ergs s~ 1 cm -2 arcsec -2 . The condensed phase of the IGM (T < 10 4 ' 5 and p/p;, > 10 3 ) is represented by the dot-dashed lines. The 
brightest He II A1640 blobs occupy a specific range of temperature (10 s < T < 10 6 ) and density (10 2 < p/pb < 10 3 ), while the brightest 
H I A1216 blobs tend to come from the high density gas. Thus the H I A1216 cooling map is greatly affected by the condensed phase of the 
IGM. To be conservative, we thus exclude the condensed phase gas particles in generating the cooling maps (Fig. |3J. 



He II A1640 luminosities as weighting factors. Thus each 
diagram represents the phases that we could actually 
probe by observing each line. In each phase diagram, 
the condensed phase of the IGM is delineated by dot- 
dashed lines. The sharp edge of the condensed phase at 
T ~ 10 4 K arises from the lack of metal-line cooling in 
our simulations. The diamonds in Figure 0] indicate the 
lines of sight that have He II A 1640 orHlA1216 surface 
brightnesses larger than 10 -19 ergs s _1 cm -2 arcsec -2 . 
Most He II A 1640 emission comes from a specific range of 
temperature (10 5 < T < 10 6 ) and density (10 2 < p/p b < 
10 3 ) that is remote from the condensed phase, whereas 
the brightest H I A1216 blobs have significant amounts 
of condensed phase gas. Because self-shielding becomes 
important in the condensed phase, we exclude this phase 
in calculating the cooling maps (in Fig. |2J) to produce 
our most conservative predictions. As we expected, this 
cutoff does not affect the He II cooling maps seriously, 
but does affect the H I A1216 cooling map dramatically, 
as gas particles with T < 10 4 5 K cannot contribute to 
He II collisional excitation cooling but only to H I colli- 
sional excitation cooling. Therefore, our predictions for 
the He II A1640 cooling radiation are far more robust 
than for H I A1216. 

The next factor that could modify the Lya cooling 
maps is the radiative transfer of Lya photons, which 
could alter the shapes and surface brightness profiles of 
the blobs substantially. Both the H I A1216 and He II 
A304 photons produced in the optically thick medium 
will be transported to the outer region by resonant 
scattering until the optical depth becomes smaller than 
r ~ 2/3. Lya photons escape eventually by scattering 
into the optically thin damping wing in the frequency 
domain, unless they are extinguished by dust. We ex- 



pect that the IGM at z ~ 3 contains little dust and 
that the cooling emission from the IGM is sufficiently far 
from the star-forming regions since we exclude the high 
density gas particles in our condensed phase cut case. 
Therefore, the net effect of the resonant scattering in the 
spatial and frequency domains is to smooth the surface 
bright ness out to the la st scattering surfaces. For exam- 
ple, iFardal et all {2001) resorted to resonant scattering 
to explain the observed size of the Steidel blobs. How- 
ever, owing to the complicated structure of the density 
and to turbulent velocity fields, it is difficult to predict 
how much radiative transfer blurs the surface brightness 
of the cooling blob. The large bulk motions will espe- 
cially affect the transfer of Lya photons. Depending on 
the optical depth and velocity field, photons can often 
undergo very little spatial diffusion and just random walk 
in velocity space unt il they reach a frequency where the 
optical depth is ~ 1 i|Zheng fc Miralda-Escudell2002j) . A 
Monte Carlo Lya radiative transfer calculation would be 
an ideal tool to make more real istic spatial and frequency 
maps of Lya cooling radiatio n ijZheng fc Miralda-Escudel 
120021 iKollmeier et al.ll2005|) . In Figure we show the 
profiles of density, temperature, velocity, and H I and 
He II emissivity for a line of sight toward a typical Lya 
blob to illustrate the complicated structure of these quan- 
tities. Because we do not take into account these radia- 
tive transfer effects in our cooling maps, the H I A1216 
cooling map (the upper middle panel in Figure should 
be smoothed to better represent reality. In contrast, be- 
cause most He II resides in the ground state, making the 
IGM optically thin to the He II A1640 line, our He II 
cooling maps should be accurate. 

Besides gravitational cooling and photoionization 
caused by the UV background, supernova feedback is an- 
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Fig. 5. — Profiles of density, temperature, H I A1216 (bold) 
and He II A1640 (light) cooling rates, and velocity for a line of 
sight. Owing to the complicated structure of the density and the 
turbulent velocity fields, it is difficult to predict how much the 
radiative transfer blurs the surface brightness of Lya cooling blobs. 
However, because the IGM is optically thin to He II, our He II 
cooling maps (Fig. [3 should be accurate. 



other energy source included in our simulations. When 
stars form in the simulations, supernova feedback energy 
is deposited into the surrounding gas in the form of heat. 
Thermal energy deposited into dense, rapidly cooling gas 
is quickly radiated away, so feedback contributes some- 
what to the cooling emission. However, we find that 
our density-temperature cutoff for the condensed phase 
effectively removes all the star-forming gas particles at 
each time step. Because the supernova thermal input 
is directly proportional to the star formation rate, our 
cooling maps without the condensed phase should not 
be seriously contaminated by supernova feedback energy. 
Thus the cooling maps in Figure |3] with the condensed 
phase removed are still robust lower limits of the flux 
from the gravitational cooling. Fardal ct al. ( 2001) also 
show that while the re-radiated supernova energy domi- 
nates at lower luminosity, gravitational cooling becomes 
the dominant source as the mass and luminosity increase. 

3.2. Properties of Cooling Sources 

To study the properties of individual H I A1216 and 
He II A1640 sources like the ones shown in the cooling 
maps, we identify discrete groups of gas particles asso- 
ciated with individual dark matter halos and then cal- 
culate the total Lya and He II A1640 luminosities for 
these sources. To find the dark matter halos, we ap- 
ply a friends-of-friends algorithm with a linking length 
that is 0.25 times the mean inter-particle separation. 
We count a gas particle as a member of the source, i.e. 
blob, associated with the dark matter halo if the dis- 
tance from the potential center is less than the virial 



radius of the halo. We then add the Lya and He II 
A1640 luminosities of the particles to obtain the total 
cooling luminosities of the blob. We restrict our analy- 
sis to blobs with more than 64 gas particles and 64 dark 
matter particles to mitigate numerical resolution effects. 
Thus the smallest halo in the 22 Mpc simulation has a 
gas mass of M gas = 6.3 x 10 9 M Q and dark matter mass 
of Mdark = 5.0 x 10 10 M Q . These masses decrease to 
M gas = 8.5 x 10 8 M Q and M dark = 6.8 x 10 9 M Q in the 
11 Mpc simulation. 

H I A1216 and He II A1640 cooling luminosities show 
tight correlations with halo mass and star formation rate 
(Fig. HJ. The open squares, crosses and circles represent 
the luminosities for the three different emissivities dis- 
cussed in H2.2I the optically thin case, the self-shielding 
correction case, and the condensed phase cut cases, re- 
spectively The correlations are as one would expect: the 
more massive a galaxy is, the more gas accretes onto the 
galaxy, resulting in more cooling radiation and a higher 
star formation rate. The distribution of cooling lumi- 
nosity is continuous, and we do not find any evidence 
that extended Lya or He II emission originates only from 
high-mass systems or high density regions. 

The solid line in the upper panel in Figure indi- 
cates the Lya emission due to recombination from stel- 
lar ionizing photons. We assume a conversion factor of 
f Lya = 2.44 x 10 42 ergs s _1 for a 1 M Q yr _1 star forma- 
tion rate with no dust absorption, no escaping ionizing 
photons, a Salpeter IMF, and solar metallicity. Under 
these assumptions the Lya emission from star formation 
in galaxies always dominates the Lya emission from the 
surrounding IGM, even in the (most optimistic) optically 
thin case. This resu lt is consistent with th e predi ctions of 
iFardal et all l|2001fl and iFurlanetto et alJ (|2005l) . 2 The 
strong correlation between the Lya cooling rate in the 
optically thin case (squares) and the star formation rate 
results from the fact that the gas in the condensed phase 
tends to satisfy the star formation criteria of the simu- 
lation and is likely to form stars in the next time step. 
In contrast, the He II emission caused by star formation 
is quite uncertain because only extremely low metallic- 
ity (Z < 10 -5 ) stars can emit the hard ionizing photons 
necessary to ionize He II. However, star formation at 
z = 2 — 3 is unlikely to be dominated by Population III 
or extremely low metallicity stars. Unlike for Lya, the 
contribution of star formation to He II must be negligi- 
ble. We discuss this point in more detail in [Q] 

Figure shows the Lya and He II luminosity func- 
tions (LFs) for the three emissivity cases at z — 2 and 3. 
The LFs include emission only from the IGM. The solid, 
dashed, and dot-dashed lines represent the optically thin, 
the self-shielding correction, and the condensed phase 
cut cases, respectively. The horizontal dotted lines in- 
dicate the number density of one halo in the simulation 
volume. Note that the distributions extend to brighter 
blobs as we increase the simulation volume because the 
larger simulations contain higher mass systems. Note 
also that the Lya luminosity function is very sensitive to 
the assumed emissivities, whereas the He II A1640 cool- 

2 In contrast, we do not find the trend of Fardal ct al. (2001) 
in which Lya from cooling radiation dominates the Lyo from star 
formation in more massive systems. We suspect that this difference 
is a consequence of including a photoionizing background in the 
simulation analyzed here. 
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Fig. 6. — Lycr and He II luminosity as a function of the halo viral mass and the star formation rate in the 22 Mpc simulation at z = 2. 
The open squares, crosses and circles represent the three different emissivity predictions: the optically thin case, the self-shielding correction 
case, and the condensed phase cut case, respectively. In the right panels, we plot only blobs with a baryonic (star + gas) mass larger than 
200 mgpjj. Below this mass limit, the derived star formation rates are not reliable owing to our limited resolution. As discussed in the 
text, Lyo luminosity changes dramatically depending on the prescription used for the self-shielded phase. The correlations between the 
cooling luminosity, halo mass, and SFR are as one generally expects: the more massive a galaxy is, the more gas accretes onto the galaxy, 
resulting in more cooling radiation and a higher star formation rate. The solid line in the upper right panel represents the Lya emission 
expected from star formation, assuming a conversion factor, /L ya = 2.44 X 10 42 ergs s~ x , for a IMq yr — 1 star formation rate with no dust 
absorption, no escaping ionizing photons, a Salpeter IMF, and solar metallicity. Note that under these assumptions the Lyo emission from 
star formation always dominates the cooling emission from the surrounding IGM. 



ing luminosity depends much less on the treatment of 
the self-shielded or condensed phase of the IGM. This 
large variatio n of the Lyq luminosity is consistent with 
the results of lFurlanetto et all lj2005fl . 

3.3. Detectability and Observational Strategy 

To estimate the detectability of cooling emission from 
the extended sources, we convert the rest-frame cooling 
maps at z — 2 and 3 into observed surface brightness 
maps and rebin them with a pixel scale of 0'.'5 x 0'.'5 to 
mimic the independent resolution elements of ground- 
based observations (Fig. 0J. Figure |H1 shows the surface 
brightness distributions of the rebinned cooling maps at 
z = 2 and 3, assuming the conservative condensed phase 
cut case. Note that the distributions depend strongly on 
the size of the bins in the surface brightness maps, be- 
cause the bright, small-scale structures are smoothed out 
by binning. We express each surface brightness distribu- 
tion in terms of the number of binned pixels per comoving 
volume and also the number of pixels per square arcmin 
if one were to observe through a R — 100 narrow band 



filter. The projected angular extents of the 11 Mpc sim- 
ulation at z — 2 and 3 are 11.4' and 9.4', respectively. 
The depth of the 11 Mpc simulation is Az ~ 0.013 and 
0.019 for z — 2 and 3, respectively. 

Deep, wide- field (~ 30' x 30'), narrow-band (R > 
100) imaging is an effective way to detect cooling 
radiation, because sky noise dominates in this low 
surface brightness range. For example, the average 
sky background at 6500A on the ground is ~ 10~ 17 

ergs s _1 cm~ 2 arcsec -2 A , comparable to our esti- 
mates for the brightest blobs. In Figure |H1 we show the 
5 a detection limits for typical R = 100 narrow band 
imaging with an 8m-class telescope and for R = 1000 
imaging with a hypothetical 30m telescope. We assume 
a peak system throughput of ~ 35%, a Mauna Kea sky 
background (for the 50% dark condition), and a 30-hour 
exposure time. We estimate the signal-to-noise ratios for 
one binned pixel (0'.'5 x O^) , which corresponds to > 2x2 
instrumental pixels in ground-based CCD detectors. 



He II Cooling Lines 



9 



V 

a 
S 



A 
C 



a 

0) 
T3 

U 
0) 

B 

3 
C 




log L (erg s ') 

Fig. 7. — Luminosity functions of Lya and He II A1640 cooling radiation at z = 2 (upper) and 3 (lower). The solid, dashed, and 
dot-dashed lines represent the optically thin, the self-shielding correction, and the condensed phase cut cases, respectively. In the upper 
panels, for each emissivity case, the LFs from the 11 and 22 Mpc simulations are denoted with light and bold lines, respectively. The 
horizontal dotted lines indicate the number density of one halo in the 11 Mpc and 22 Mpc simulations. Note that the 22 Mpc results 
extend to higher luminosity than the 11 Mpc results, and that even the larger simulation may underestimate the density of high luminosity 
systems. In the lower panel s, we show onl y results from the 22 Mpc simulation. For comparison, we show the LF of the Subaru Lya blob 
sample (the crosses; Matsuda ct al. 2004). The Subaru Lya blobs could arise from a variety of mechanisms (see text), including cooling 
radiation. It is suggestive, however, that their LF lies between our predictions for the self-shielding correction and condensed phase cut 
cases. Note also that the Lya LF is very sensitive to the assumed emissivities, whereas the He II A1640 cooling luminosity depends much 
less on the treatment of the self-shielded or condensed phase of the IGM. 



3.3.1. H I A1216 

We predict that H I A1216 cooling emission from 
the brightest blobs at z — 2 and 3 is detectable by 
6-8m class telescopes with moderate resolving power 
(R = 100). The limiting sensitivity of current surveys 
for high-z Lya emitters is ~ 10~ 18 ergs s _1 cm~ 2 (e.g., 
Malh otra &: Rhoadsl 12004 and references therein) . It 
is encouraging that even our most conservative predic- 
tions suggest that the Lya blobs arising from gravita- 
tional cooling radiation are detectable with a reasonable 
amount of telescope time. 

The Lya surface brightness of the largest system in 
our z = 3 simulation (M ha io ~ 6.5 x 1O 12 /i -1 M ), corre- 
sponding to the brightest blob in Figure |3 is consistent 
wit h the mean surface b rightnesses of the Lya blobs of 
the lMatsuda et all l)2004[) sample (represented with small 
vertical bars in the Figure |8J . Note that our predicted 
Lya blob luminosities depend on the different emissivi- 
ties for the self-shielded phase and/or the exact location 
of our density-temperature cut of the condensed phase. 
Though the surface brightnesses of the predicted and ob- 
served blobs are consistent, the luminosity of our bright- 
est Lya blob is fainter than that of observed blobs (see 
Fig. UJ, possibly because of our conservative assumptions 
for the self-shielded phase. 



3.3.2. He II A1640 

A pixel-by-pixel comparison of the H I A1216 and He II 
A1640 cooling maps reveals that, without the condensed 
phase, the He II A1640 flux is always > 10 x fainter than 
that of H I A1216. The He II A1640 emission could be 
1000 x fainter than Lya in the optically thin case, i.e., 
the most optimistic Lya prediction. Nonetheless, detec- 
tion of the He II A1640 cooling line from z = 2 sources is 
clearly feasible with 6-8 meter class telescopes, and even 
possible at z = 3. Though the number statistics of bright 
blobs are limited by the relatively small volume of our 
simulations, we expect one source in the 11 Mpc simula- 
tion and six sources in the 22 Mpc simulation at z = 2 
with areas of > 0'.'5 x 0'.'5 above the surface brightness de- 
tection threshold of ~ 5 x 10~ 18 ergs s _1 cm -2 arcsec -2 
(R = 100 arrow in Figure 8). Thus the space den- 
sity of the sources from which we could detect not only 
Lya but also He II A 1640 emission with narrow band 
imaging corresponds to a comoving number density of 
- 5-7x 10~ 4 h 3 Mpc~ 3 or - 0.02 arcmin" 2 per R = 100 
filter (Az = 0.03; AA ~ 49A for He II A1640). 

If we consider the larger survey volume typically acces- 
sible by modern narrow band imagers, we could expect 
better survey efficiency than that described above. Be- 
cause cosmological simulations do not contain power on 
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scales larger than their finite sizes, the largest objects 
in the simulation are typically underestimated in num- 
ber and size. For example, the luminosity functions and 
surface brightness distributions in Figs. and |3 extend 
their bright limits as the volume of the simulation in- 
creases. Therefore, one might expect that the detection 
of even brighter blobs by surveying more volume. For 
example, many current wide-field imagers and spectro- 
graphs have half degree field of views, so it is possible to 
survey a volume ~ 17 times larger than encompassed by 
our 11 Mpc simulation at z — 2 (or a volume ~ twice 
that of the 22 Mpc simulation). Therefore, if we naively 
extrapolate the number of detectable He II sources in our 
11 Mpc and 22 Mpc simulations, then 17 (±17) and 13 
(±5) He II sources would be detected, respectively. The 
numbers within parentheses indicate Poisson errors. 

Our results predict that bright He II sources are always 
bright Lya cooling blobs. Observationally, the difficul- 
ties in searching for He II cooling sources could be eased 
(1) by pursuing narrow-band Lya imaging first, detect- 
ing Lya blobs, and looking for He II A1640 emission in 
those blobs with follow-up observations, or (2) by adopt- 
ing a combine d multislit spectroscopy + narrow-band fil- 
ter approach l|Martin & Sawickil200l iTran et alJl2004^ 
to identify Lya blobs (which can then be targeted for 
Hell). 

The latter technique is potentially quite effective de- 
spite the faint surface brightness of the Lya and He II 
blobs. This technique employs multiple parallel long slits 
with a narrow band filter that limits the observed spec- 
tral range to a few hundred angstroms and, by dispers- 
ing the sky background, achieves better sensitivity than 
narrow-band imaging alone. In the sense that this tech- 
nique trades off survey volume (or sky coverage) to go 
deeper in flux, it is about as efficient as simple narrow 
band imaging for surveys of Lya emitters, which are not 
as extended as Lya blobs. However, the multi-slit win- 
dow technique is superior to narrow band imaging for 
low surface brightness objects like the Lya and He II 
blobs discussed here. It has the further advantages that 
1) it provides the spectral and kinematic data necessary 
to distinguish the origins of blob emission (<2J) and 2) it 
enables us to exclude contaminating emission lines, such 
as Ha, H/3, [O III], and [O II], from nearby star- forming 
galaxies by measuring the line shapes (e.g. line asymme- 
tries and line doublets) and blueward continuum. If this 
technique is employed with large field-of-view imagers, 
the survey volume covered by the multiple slits is still 
reasonably large (~ 10% of the whole field of view). 

It is possible to search for He II cooling radiation at 
lower redshifts than z ~ 2 — 3. For example, z rs 1.5 
is the lowest redshift at which He II A 1640 still lies at 
an optical wavelength (A f, s ~ 4100 A). Because metal 
abundanc es in the IGM do not change very much over 
z ~ 2 — 4 ( Sch ave et al.l20 03i) . it is unlikely that our basic 
assumption of a primordial composition (j}2J is violated 
seriously at z sw 1.5. Any blind search for He II A 1640 
blobs at z = 1.5 — 2 will be contaminated by H I Lya 
emission from z « 3 sources if only one emission line is 
identified in the spectrum. In this case, the blob's red- 
shift could be further constrained by obtaining a redshift 
for the galaxy it surrounds. 

He II A1640 cooling radiation at very low redshifts 



(z < 0.5) is potentially detectable in t he ultraviolet us- 
ing UV satellites or HST. For example, Furlanetto et 
(2003) show that detection of the bright cores of H I 
A1216 emission from z < 0.5 sources is feasible with deep 
wide-field UV imaging, e.g. with The Galaxy Evolution 
Explorer (GALEX) or the proposed Space U ltraviolet- 
Visible Observatory fSUVO: IShull et al] li999L Because 
H I A1216 and He II A1640 trace different phases of 
the shown in Figure 0] combined observations 

of these two lines, e.g., of their morphologies and line 
ratios, would probe different phases of the IGM. 

3.3.3. He II A304 

In contrast to H I A1216 and He II A1640, He II A304 
photons redshifted to wavelengths shorter than the pho- 
toionization edge of H I and He I (912A and 504A, respec- 
tively) can be absorbed by neutral hydrogen and neutral 
helium. Even if they escape the blobs, He II A304 pho- 
tons are removed from the line of sight owing to cumula- 
tive absorption by the intervening neutral IGM, includ- 
ing the Lya forest and damped Lya systems. We esti- 
mate the transmission of He II A304 through the inter- 
ve ning IGM using Monte Ca rlo simulations as described 
in lMoller fc JakobserJ l)1990|) with the updated statistics 
of Lyman forest and Lyman limit systems (i.e., num- 
ber density evo l ution and column density distribution) 
from pakobscn] l|1998fl . For the emitters at z = 3, we 
find that the transmission factor averaged over all lines 
of sight is - 12% and that 67% (76%) of the sightlines 
will have transmission lower than 1% (10%). Therefore, 
though the He II A304 emissivity is roughly 10 x higher 
than that of He II A1640 (Fig. [l]in ffH|) . we expect the 
He II A304 cooling map to be fainter than that of He II 
A 1640 in most cases and to vary strongly from sightline 
to sightline. 3 

He II A304 photons can also be destroyed by the H I 
and He I inside the blobs, because He II A304 photons 
experience a large number of scatterings before escaping. 
The destruction probability by H I and He I atoms per 
scattering is given by 

rtHiCHi + "Hci^Hoi , n s 

e = , (3) 

"Hi^Hi + TlHeiOHel + "HellOLya 

where cthi and anei are the photoionization cross sec- 
tions of H I and He I at 304 A, respectively, and CLya is 
the integrated scattering cross section of He II A304. The 
abundances of H I and He I atoms are smaller than for 
He II, and their photoionization cross sections are also 
much smaller than the resonant cross section of He II 
A304 by a factor of < 2 x 10~ 5 . We estimate that the 
destruction probability is e > 5 x 10 -8 at a tempera- 
ture of T ~ 10 48 without applying the self-shielding cor- 
rection. In the self-shielded regions where more neutral 
hydrogen can reside, this probability rises. Thus the es- 
cape probability of a He II A304 photon from a blob is 
/igm — (1 — c) Nt , where N T ~ t 2 is the number of 
scatterings required to escape the blob, if it is approxi- 
mated by an optically thick slab. For example, we obtain 
/igm — 0.7% for THen = 10 4 . Therefore, we cannot ig- 
nore the absorption by H I and He I atoms in the high 

3 If one sightline does not have any Lyman limit or damped Lyo 
systems, it will have ~87% transmission on average due only to 
Lyman forest systems. 
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Fig. 8. — Distributions of surface brightnesses in the cooling maps for the IGM at z = 2 and 3. Bold and light solid lines represent the 
surface brightness histograms of rebinned (0"5 X 0'.'5) pixels for He II A1640 and H I A1216, respectively, in the 11 Mpc simulation. We 
show the S/N > 5 detection limits for a 30- hour observation with an 8- meter telescope and R = 100 narrow-band filter, and with a 30-meter 
telescope and R = 1000 filter. The right y-axis represents the number of binned pixels per square arcminute per the redshift width for the 
R = 100 filter. The dotted lines at the bottom of each panel indicate one detection in the simulation. The dot-dashed lines represent the 
surface brightness distributions for the 22 Mpc simulation. Detection of the brightest H I A1216 and He II A1640 blobs at z = 2 and 3 
is possible with deep narrow-band imaging on a 6-8 meter class telescope. Note that the larger simulation extends the bright tail of the 
surface brightness distribution. Even larger simulation volumes might therefore predict still brighter, thus more easily detectable, systems. 



density gas. However, the number of scatterings iV r is 
very difficult to estimate correctly because of the com- 
plex density and velocity structure of the blob, unless one 
carries out full 3-D hydro-radiative transfer calculations 
(which are beyond the scope of this paper). For certain 
geometries and velocity fields, bulk motions of the gas 
will help the He II A3 04 photons escape the blob wi th 
fewer scatterings (e.g., Zheng & Miralda-Escude 2002). 

Owing to intervening absorption and the destruction 
inside the blobs, He II A304 is the most uncertain cool- 
ing line we consider. Although He II A304 is dimin- 
ished significantly by the intervening IGM, if the es- 
cape fraction from the IGM is significant (/igm — 1), 
the detection of He II A304 may not be out of ques- 
tion with a large apert ure UV /optical op timized space 
telescope (e.g., SUVO: IShull et al.lli"999|) . An advan- 
tage of observing He II A304 in the far ultraviolet in 
space is that the sky background is very low (~ 10~ 23 
ergs s _1 cm~ 2 arcsec -2 A -1 at 1250 A) compared to the 
optical (~ 10~ 18 ergs s _1 cm~ 2 arcsec~ 2 A -1 at 6500 
A), except for the geocoronal emission lines (e.g., Lya 
1216 A and O I 1304 A). Thus, if these geocoronal emis- 
sion lines could be eliminated with blocking filters or by 
adopting an L2 orbit, the direct detection of He II A304 
is feasible. In this case, the detector noise — especially 
the dark current — will dominate. Recent developments 
in UV detector technology are very promising, so the 
possibility of studying these blobs at those wavelengths 
remains open. 

4. DISCUSSION 



Until now, we have only considered the cooling ra- 
diation from gas that is losing its gravitational energy, 
falling into a galaxy-sized dark halo, and ultimately form- 
ing stars. Photoionization by these stars is another pos- 
sible heating source for the blobs. Starburst-driven su- 
perwinds or AGNs, which are not included in our sim- 
ulations, are other p otential blob energy sources (e.g., 
see t he discussions in iSteidel et aTll2000l iMatsuda et all 
2004). Although the radiation from gas heated by these 
feedback processes is not generally termed "cooling ra- 
diation", the energy injected into the surrounding gas 
can also be released through line emission. Thus our 
estimates for cooling emission might be lower limits for 
the actual fluxes in the cooling lines. In this section, 
we assess whether other H I A1216 and He II A 1640 
sources overwhelm our gravitational cooling signals and 
then discuss how to discriminate among these other pos- 
sible mechanisms in order to use H I A1216 and He II 
A1640 cooling lines to study gas infall into galaxies. 

4.1. Photoionization by Stellar Populations 

UV photons from massive stars in a galaxy or blob 
ionize the surrounding interstellar medium, and the re- 
combination lines from these nebulae could contribute to 
the Lya and He II line fluxes. Generally, the recombina- 
tion line luminosity is proportional to the star formation 
rate (SFR) and is given by 

(SFR \ 
M yr" 1 j 

where Td us t is the dust optical depth for the ionizing 
continuum in the interstellar medium (ISM), / osc is the 
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fraction of ionizing photons that escape the star-forming 
galaxy, /igm is the fraction of photons that escape the 
surrounding IGM, and /n nc is the conversion factor from 
the SFR to the line luminosity in ergs s^ 1 . This SFR 
conversion factor depends on the metallicity, initial mass 
function (IMF), and evolutionary history of the stars in 
the blobs (e.g., a lower metallicity and a top heavy IMF 
produce more ionizing photons and thus more recombi- 
nation line photons). 4 

The conversion factor for H I A1216, /1216, is large 
enough to make it difficult to distinguish the cooling 
lines (of IGM origin) from the r ecombination-in duced 
lines (of ISM origin). For example. ISchaereil l)2003l) finds 
/1216 = 2.44 x 10 42 ergs s _1 for a constant star formation 
history, solar metallicity, and a Salpeter IMF with a mass 
range of 1 - 100 M s . Thus for a SFR = 10 M Q yr" 1 , 
/ csc ~ 0.1, e -Tdust ~ 0.1, and /igm — 1, we obtain the 
observed flux, F1216 — 2.9 x 10 -17 ergs cm" 2 s" 1 , due to 
the stars in a Lya blob at z = 3, which is comparable 
to the surface brightness of the brightest blobs in our 
simulations (see fc>.3l for the blob luminosity functions). 
Therefore, the contamination of the H I A1216 line by 
stars is not negligible, unless the Lya photons from star- 
forming regions are heavily absorbed by a dusty ISM 
(e.g., in highly obscured sub- millimeter galaxies or Ly- 
man break galaxies with the damped Lya absorption). 
Because Lya cooling radiation is produced sufficiently 
far from the star-forming region and thus should be less 
susceptible to dust attenuation than the Lya emission 
from the stellar populations, it might be possible to iso- 
late extended Lya cooling radiation in these galaxies. 
However, in the case that Lya photons emitted by stars 
escape the galaxy (e.g., Lya emitters), it will be chal- 
lenging to distinguish the Lya cooling radiation from the 
Lya produced by the stellar populations unless the var- 
ious parameters such as / eS c, /igm, and SFR are fully 
constrained. 

In contrast, He II A1640 emission appears to be lim- 
ited to very small metallicities (log(Z/Z ) < —5.3) and 
Population III objects, because stars of solar or sub- 
solar metallicities emit few if any He I I ionizing pho- 
tons (iBromm. Kudritzki. fc Loebl 120011 iSchaererl Hool 
iTimlinson. Shull. & Venkatesanl 120031) . Using He II 
A1640 to detect the first hard-ionizing sources such as 
metal- free stellar populations, the first miniquasars, or 
even stellar populatio ns before the reionization epoch has 
been proposed (e.g., J^mfingoji^G.ir^ux." fc Shnll f2001: 
lOh. Haiman. fc ReeslEoOlt iBarton et al.ll2004D . In this 
paper, we take advantage of this fact to discount the 
contributions of stellar populations to the He II A1640 
cooling line. Even for Z = 10~ 5 , /i640 = 1-82 x 
10 40 ergs s" 1 (~ 6 x 10~ 4 /i2i6) for an extremely top 
heavy IMF containing only stars in the range 50 — 
500 Mq. For the same assumptions used in the Lya 
calculation above, we obtain an observed flux of -F1640 — 

4 Note that the equation does not take into account the 
dust absorption of Lya photons by ISM after Lya photons escape 
from the H II region. When the H I optical depth in the ISM is 
high enough, Lya photon will experience a large number of res- 
onant scatterings before escaping the galaxy. These scatterings 
will increase the effective dust optical depth and destroy Lya pho- 
tons preferentially. However, the actual optical depth is likely to 
strongly depend on the kinematics of neutral gas and the geometry 
of the galaxy, which are hard to quantify. 



2.2 x 10 -19 ergs cm -2 s _1 from the blob stars at z = 3, an 
order of magnitude below the brightest He II A1640 blobs 
in our simulations. Therefore, it is very unlikely that 
the He II A1640 photons originating from stars contami- 
nate the gravitational cooling emission, unless significant 
numbers of metal-free stellar populations are forming. 
Thus He II A1640 cooling radiation is much less contami- 
nated than H I A1216 by recombination lines originating 
from star-forming galaxies. 

The only caveat is the possibility of He II A1640 emis- 
sion arising not from stars directly, but from the hot, 
dense stellar winds of Wolf-Rayet (W-R) stars, the de- 
scendants of O stars with masses of M > 20 — 30 Mq. 
W-R populations formed in an instantaneous starburst 
at high redshifts would not seriously contaminate the 
He II A1640 cooling radiation, because W-R stars are 
very short-lived (< 3 Myr) and their number relative to O 
stars (W-R/O) drops as the metallicity decreases below 
solar. In the case of continuous star for mation, a stellar 
popu lation synthesis model (Starburst99; Lei therer et alJ 
1999) predicts that the maximum number of W-R stars is 
reached ~10 Myr after the initi al burst. Using the He II 
A 1640 luminosity of a W-R star l|Schaerer fc Vaccall99S() 
and the number evolution of W-R stars from Starburst99 
under the assumptions of a Salpeter IMF (1 — 100 Mq), 
sub-solar metallicity [Z ^ 0AZ@), and a massive SFR of 
100 Mq yr _1 over at least 10 Myr, we estimate the He II 
A1640 line luminosity due to W-R stellar winds to be 
< 10 42 ergs s _1 . Although this He II A1640 luminosity 
is comparable to the predicted He II cooling radiation, it 
is possible to discriminate between the two He II A1640 
sources in individual objects because the emission from 
W-R winds should be much broader (e.g. ~ 1000 km s _1 ; 
see Fig. [5] in ^4.31 which is relevant here even though it 
is presented in the context of the galactic superwind sce- 
nario) . 

One way to test our predictions in this section is to 
look more closely at the He II A1640 emission associ- 
ated with high rcdshift star-forming galaxies, i.e., the Ly- 
man break galaxies (LBGs) with vigorous star formation 
rates. IShaplev et alJ ((2003) show that composite spectra 
of LBGs have very broad (FWHM - 1500 km s" 1 ) He II 
A1640 profiles regardless of their Lya emission strength. 
While they attribute the He II features to W-R stel- 
lar winds, those authors have difficulty reproducing the 
strength of the He II lines using stellar population syn- 
thesis models with reasonable parameters. Because of 
this inconsistency, we speculate that some fraction of the 
He II features may come from the cooling of gas falling 
into these galaxies along line of sight. It would be worth- 
while to obtain high signal-to-noise spectra of individual 
LBG's and their surroundings to see if the He II line 
is present, especially outside the galaxy, and relatively 
narrow. 

4.2. Photoionization by AGNs 

AGNs inside the star-forming regions of blobs could 
photoionize the surrounding gas and generate He II 
A1640 as well as H I A1216 emission. The predicted size 
(a few arcseconds) and surface brightness (~ 10~ 18 — 
10~ 16 ergs s _1 cm~ 2 arcsec~ 2 ) of an extended Lya blob 
enshroudin g a quasar are consis tent with the observed 
quantities ijHaiman fc ReesH200lD . 

How many AGN-powered sources might we expect in 
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a Lya/He II blob survey? Unfortunately, it appears that 
there is no easy way to predict Lya/He II luminosity 
from the surrounding IGM, because we do not know 
how much neutral IGM is distributed around the AGN. 
Therefore we take a conservative approach to estimate 
the number of the AGN-powered sources. First, we es- 
tablish a simplistic relationship between the induced Lya 
(or He II) blob luminosity and the X-ray luminosity of 
the AGN, then we estimate the number of AGN-powered 
blobs based on the known hard X-ray luminosity func- 
tion of AGN at z = 2 — 3. If there are not many relative 
to the number of blobs powered by gravitational cooling 
radiation, then we could conclude that they are unlikely 
to complicate the interpretation of extended Lya/He II 
sources. 

We assume that all the ionizing photons from an AGN 
with a simple power-law spectrum, L v = L§(y / vo) a , are 
absorbed by the surrounding medium and re-emitted as 
recombination lines. The line (Lya or He II) luminosity 
of the surrounding blob is then given by: 

t n- r L °{ v Y ™ 

-Mine — ClinoW — Qine / — I , 10) 

where Q is the number of ionizing photons emitted per 
unit time, vll is the frequency of the Lyman limits for 
the hydrogen and He II \hv LL — 13.6 eV and 54.4 eV, 
respectively), and the line emission coefficient cij nc in 
ergs is the energy of the line photon emitted for each 
H I or He II ionizing photon. For case-B recombination 
with an electron temperature of T e = 30, 000K and an 
electron number density of n e — 100 cm~ 3 , we obtain 
Ci2ifi = 1-04 x 10" 11 and c 16i0 = 5.67 x 10~ 12 ergs (c.f. 
ISchaererll2003|) . We adopt a spectr al index a = —1.8 for 
the extreme UV l|Telfer et al.l2002|) and assume that this 
a is valid even in the X-ray. The hard X-ray luminosity 
of AGN (Lx) is simply given by the integration of L v 
between 2 keV and 8 keV. 

Once the Lya (or He II) luminosity is monotoni- 
cally linked with the AGN X-ray luminosity, we esti- 
mate the number density of AGN-powered sources from 
the h ard X-ray lu minosity function (e.g., iBarger et al.1 
120031 iCowie et al.1l2003|) . To power a blob with L Lya 
= 10 43 ergs s" 1 , an AGN must have L x > 10 418 
ergs s" 1 , which would generate a He II blob with Lhch > 
10 41 - 7 ergs s" 1 . A round the required X-ray luminosity, 
ICowie et al.1 (|2003) derive the number density of X-ray 
selected AGNs regardless of their optical AGN signa- 
tures to be 1.3 x 10~ 5 < <Z>(L X > 10 42 ergs s" 1 ) < 
1.4 x 10~ 4 Mpc~ 3 at 2 < z < 4. The extreme upper limit 
was determined by assigning all the unidentified sources 
in the survey to 2 < z < 4 and is thus very conserva- 
tive. For the brightest cooling sources in our simulations 
(Fig. CJ, we find $(L LyQ > 10 43 ) ~3x 10~ 5 Mpc" 3 and 
~ 9 x 10 -4 Mpc~ 3 for the condensed phase cut and the 
self-shielding correction cases, respectively. 5 Note again 
that our assumption that all the ionizing photons from 
all AGNs are absorbed to produce Lya/He II photons is 
very conservative and that we are clearly over-predicting 
the number of AGN-powered blobs. However, even un- 
der this conservative assumption the number density of 

5 The number de nsity of Lya blobs (of unknown origin(s)) from 
the Subaru survey ( Matsuda et al. 2004) lies between these two 
cases. 
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Fig. 9. — Distribution of He II A 1640 flux- weighted velocity dis- 
persion of the gas particles associated with individual dark matter 
halos. The different lines represent the velocity dispersions in the 
x, y, and z directions. Note that the velocity dispersion due to 
gas accretion is less than 400 km s , in contrast to the typica l 
galactic wind speed of ~ 400-800 km s — 1 iHeckman et al.ll200fl1 . 
Thus the width of the optically thin He II A1640 line could be used 
as a diagnostic to discriminate between the galactic wind and the 
gravitational cooling hypotheses for powering Lya blobs. 



AGN-powered sources is only marginally comparable to 
the number density of cooling sources in our simulations. 
Therefore, at present, we conclude simply that a survey 
for extended Lya and He II cooling radiation is not likely 
to be swamped by AGN-powered sources. 

The above arguments are statistical, whereas distin- 
guishing gravitational cooling radiation from the emis- 
sion of AGN-photoionized gas for an individual source 
requires a multi-wavelength approach. It is therefore use- 
ful to target fields with a large amount of ancillary data 
(e.g., deep broad-band or X-ray imaging) to make an un- 
ambiguous detection of a true He II cooling blob. First, 
searching for the C IV (1549 A) emission line in the op- 
tical spectrum of the source is a good way to identify 
an AGN (e.g.. lKeel et al]ll999[) . The composite spectra 
of optically selected qu asars show bright C IV lines, but 
much fainter He II lines ijTelfer et al.12 002'). The recently 
discovered Lya blob associated with a luminous mid- 
infrared source l)Dev et al. 2005) shows unusually strong 
He II A1640 lines and C IV lines in a localized region 
near the center of nebula, suggesting that this Lya blob 
is powered, at least in part, by an obscured AGN. On 
the other hand, the absence of a C IV line in a spectrum 
with strong He II emission might indicate gravitational 
cooling gas like that in our simulations. As we discuss in 
the next section, the kinematics of the He II line can fur- 
ther constrain the origin of the He II emission. Second, 
if the AGN is heavily obscured, deep X-ray imaging of ~ 
Msec will provide the most direct probe, because X-rays 
from the AGN can penetrate the large column densities 
of gas and dust. 

4.3. Superwinds 

Alternatively, iTaniguchi fc Shioval {2000) suggest that 
galactic superwinds driven by starbursts could power the 
extended Lya blobs. In this scenario, the collective ki- 
netic energy of multiple supernovae is deposited into the 
surrounding gas, producing a super-bubble filled with hot 
and high-pressure gas. If the mechanical energy over- 
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comes the gravitational potential of the galaxies, this 
metal-enriched gas blows out into the surrounding pri- 
mordial IGM and evolves into superwinds. 

While the luminosity and sizes of the observed blobs 
are roughly consistent with the predictions of simple 
wind models, the mechanism to convert the mechanical 
energy into Lya emis s ion is not clear. Using a fast-shock 
model, iFrancis et alJ l)2001|) show that if the shocks are 
radiative, the emission from the excited gas in the shock 
itself and the photoionized precursor region in front of the 
shock can explain the observed Lya surface brightness of 
the blobs. For example, if we adopt the fiducial model of 
the pre-run sho ck grids fromlAllen et alJ J2004D (MAP- 
PINGS code bv lDopita fc Sutherland! 1996(1 with a shock 
velocity of 700 km s _1 , a number density of 10" 2 cm" 3 , 
a magnetic parameter of B/y/n = 2/xG cm 3 / 2 , and so- 
lar metallicity, then the Lya and He II A1640 emissivity 
from the shock + precursor region will be ~ 0.04 and 
0.002 ergs s" 1 cm" 2 , respectively. If the shock is per- 
pendicular to our line of sight, we would expect sur- 
face brightnesses of F^ ya ~ 3 x 10" 18 and Fn CII ~ 
1.5 x 10 -19 ergs s _1 cm -2 arcsec -2 at z = 3. This Lya 
surface brightness is roughly consistent wi th the observed 
mean surface brightness of Lya blobs (Matsuda et al. 
2004), but is an underestimate if we consider that Lya 
emission can be suppressed by various factors such as 
self-absorption and that the density of the IGM is pos- 
sibly lower than the assumed density. If the IGM in the 
preshock region is composed of neutral primordial gas, 
the UV photons produced in the post-shock plasma will 
ionize the preshock region, and the lack of an effective 
cooling mechanism other than the atomic hydrogen and 
helium lines can boost the Lya and He II A 1640 line 
emissivities significantly. However, the low metallicity 
shock grid is not currently available, and the density of 
the IGM in the preshock region and the effect of mixing 
between the metal-enriched winds and the pristine IGM 
are quite uncertain. Thus it is difficult to predict how 
much mechanical energy is released through the Lya or 
He II A1640 lines in the superwind model. 

One important feature of the superwind shock model 
is that it also predicts many UV diagnostic lines (e.g., 
C IV A1549) that have been used to study the energetics 
of the narrow-line region in AGNs. The debate about 
the origin of the Lya blobs arises mainly because Lya is 
not a good diagnostic line to discriminate between AGN 
photoionization and superwind shock-excitation owing to 
its sensitivity to resonant scattering and obscuration by 
dust. Ideally, line ratios (e.g., He II/C IV, once detected) 
could be used to discriminate among the different mech- 
anisms. 

The kinematics of the blob is potentially another test of 
the superwind hypothesis, because of the expected bipo- 
lar outflow mo t ion of the expanding shell. For exa mple, 
lOhvama et"aTl (|2003fl claim that Blob 1 of iSteidel etall 
(2000) shows both blueshifted and redshifted compo- 
nents (~ ±3000 km s" 1 ) in the central region, and they 
attribute these profiles to the expanding bipolar motion 
of a shocked shell driven by a superwind. On the other 
hand, using integral field spectrograph data. lBower et all 
l|2004() argue that Blob 1 has chaotic velocity structures 
that can be explained by the interaction of slowly rising 
buoyant material with cooling gas in the cluster poten- 



tial, and that a powerful collimated outflow alone ap- 
pears inconsistent with the lack of velocity shear across 
the blob. 

We show the He II A1640 luminosity-weighted velocity 
dispersions of the gas particles associated with individual 
blobs in Figure^] The effect of Hubble expansion and pe- 
culiar motion is included in the velocity dispersion calcu- 
lations but the thermal broadening for each gas particle is 
not. Most halos have velocity dispersions smaller than ~ 
400 km s , compared to the typical sup erwind speed of 
several hundreds to a 100 km s _1 fe.g-. lHeckman et alJ 
120001 IPettini et al1l200lH . For the superwind case, be- 
cause the observed Lya emission comes mainly from the 
shock between the outflow from a galaxy and the sur- 
rounding pristine IGM, we expect the He II emission to 
be as extended as the observed Lya emission. Therefore, 
if we observe a spatially resolved Lya and He II emit- 
ting blob, and its velocity dispersion is larger than 400 
km s , it is possible to exclude cooling radiation as the 
source of that blob. Note that there would be no ambi- 
guities in measuring the size and line broadening because 
He II A1640 is optically thin. Thus He II A1640 is a finer 
tool than H I A1216 to study the kinematic properties of 
Lya blobs. 

5. CONCLUSIONS 

In this paper, we use high resolution cosmological sim- 
ulations to study the gravitational cooling lines arising 
from gas accreted by forming galaxies. Because baryons 
must radiate thermal energy to join a galaxy and form 
stars, accreting gas produces extended H I A1216 emis- 
sion (a "Lya blob") surrounding the galaxy. We also 
expect cooling lines from singly ionized helium such as 
He II A1640 to be present within Lya blobs. We investi- 
gate whether three major atomic cooling lines, H I A1216, 
He II A1640, and He II A304 are observable in the FUV 
and optical. We discuss the best observational strate- 
gies to search for cooling sources and how to distinguish 
them from other possible mechanisms for producing Lya 
blobs. Our principal findings are: 

1. H I A1216 and He II A1640 (He II Balmer a) cooling 
emission at z = 2 — 3 are potentially detectable with 
deep narrow band imaging and/or spectroscopy from the 
ground. He II A304 will be unreachable u ntil a large 
apert ure UV space telescope (e.g. SUVO; :Sh ull et alJ 
1999) is available. 

2. While our predictions for the strength of the H I 
A1216 emission line depend strongly on how to handle 
the self-shielded gas, our predictions for the He II A1640 
line are rather robust owing to the negligible emissivity 
of He II for the self-shielded IGM below T ~ 10 4 5 K. 

3. Although He II A1640 cooling emission is fainter 
than Lya by at least a factor of 10 and, unlike Lya 
blobs, might not be resolved spatially with current ob- 
servational facilities, it is more suitable to study gas ac- 
cretion in the galaxy formation process because it is op- 
tically thin, less sensitive to the UV background, and less 
contaminated by recombination lines from star-forming 
galaxies. 

4. To use the H I A1216 and He II A1640 cooling 
lines to constrain galaxy formation models, we first need 
to exclude the other possible mechanisms for produc- 
ing Lya blobs. First, because He II A1640 emission 
from stars is limited to stars with very low metallici- 
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ties (log(Z/Z Q ) < —5.3) and Population III objects, its 
detection, unlike Lya, cannot be caused by stellar pop- 
ulations. Second, the kinematics of the He II A1640 line 
can distinguish gravitational cooling radiation from a sce- 
nario in which starburst-driven superwinds power Lya 
blobs, because the He II line width from cooling gas is 
narrower (a < 400 kms -1 ) than the typical wind speeds 
(which are factors of several higher). Third, if some frac- 
tion of the He II emitting blobs are powered by AGN, 
additional diagnostics such as the C IV line and/or X- 
ray emission can be used to discriminate gravitationally 
cooling blobs from those powered, at least in part, by 
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